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Abstract. We illustrate how recently improved low-redshift cosmological 
measurements can tighten constraints on neutrino properties. In particular we 
examine the impact of the assumed cosmological model on the constraints. We 
first consider the new HST Hq = 74.2it3.6 measurement by Riess et al. (2009) and 
the (T8(r2™/0.25)° ''i = 0.832 ± 0.033 constraint from Rozo et al. (2009) derived 
from the SDSS maxBCG Cluster Catalog. In a ACDM model and when combined 
with WMAP5 constraints, these low-redshift measurements constrain < 0.4 

eV at the 95% confidence level. This bound does not relax when allowing for 
the running of the spectral index or for primordial tensor perturbations. When 
adding also Supernovac and BAO constraints, we obtain a 95% upper limit of 
i^v < 0.3eV. We test the sensitivity of the neutrino mass constraint to the 
assumed expansion history by both allowing a dark energy equation of state 
parameter w ^ —1 and by studying a model with coupling between dark energy 
and dark matter, which allows for variation in w, fij., and dark coupling strength 
5. When combining CMB, Ha and the SDSS LRG halo power spectrum from 
Reid et al. 2009, we find that in this very general model, X]™" < O.bl eV with 
95% confidence. If we allow the number of relativistic species N^^i to vary in a 
ACDM model with J^m^ = 0, we find N^^i = ^-^^tl^lt^^X) for the 68% and 
95% confidence intervals. We also report prior-independent constraints, which are 
in excellent agreement with the Bayesian constraints. 



1. INTRODUCTION 

Atmospheric and solar neutrino experiments have demonstrated that neutrinos have 
mass, implying a lower limit on the total mass of 0.056 eV [Tj. Ongoing and future 
direct experiments will be sensitive to the absolute neutrino mass scale: for instance 
KATRIN will measure or constrain the electron neutrino mass to the 0.2 eV level 
A thermal neutrino relic component in the universe can impact both the expansion 
history and growth of structure. Cosmological probes of neutrinos are therefore 
complementary to direct experiments and provide some of the tightest constraints 
on both the number of relativistic species present in the early universe and the sum 
of neutrino masses [T]. However, up to date, these constraints were dependent on the 
assumed cosmological model. For example, bounds on the sum of the neutrino masses 
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could be relaxed by, for example, allowing coupling in the dark sector [2J or allowing 
hadronic axions [3]. In this paper we summarize how cosmic neutrinos affect cosmic 
history, examine how recently improved cosmological measurements of the low redshift 
universe can help constrain these neutrino properties, and scrutinize the impact of our 
assumptions about the cosmological model. 

In the standard model for particle physics there are three massless neutrino 
species with weak interactions. Neutrinos decouple early in cosmic history, when 
the temperature was T ^ 1 MeV, and thereafter contribute to the relativistic energy 
density with an effective number of species, Ncs — 3.046 [T]. Cosmology is sensitive 
to the physical energy density in relativistic particles in the early universe ojrei which, 
in the standard model (for cosmology), includes only photons and neutrinos io^ci = 
ujj + NoffUj^ where uJj is the energy density in photons and uj^ is the energy density in 
one active neutrino. As uj-y is extremely well constrained, Wroi can be used to constrain 
neutrino properties; deviations from Ngg = 3.046 would signal non-standard neutrino 
features or additional relativistic relics. A^eff impacts the big bang nucleosynthesis 
epoch through its effect on the expansion rate. Therefore measurements of the 
primordial abundance of light elements can constrain NeS [H O O [7] . These constraints 
rely on physics at the time of big-bang nucleosynthesis (T ~ MeV), and in several non- 
standard models the energy density in relativistic species can change at later time (e.g. 
at the last scattering, T ^ eV). Free-streaming relativistic particles affect the Cosmic 
Microwave Background (CMB) through their relativistic energy density, which alters 
the epoch of matter-radiation equality, and through their anisotropic stress [3 |9]. 
Because the redshift of matter-radiation equality is well constrained by the ratio of 
the third to first CMB peak height, the first effect defines a degeneracy between Vlmh? 
and A^eff CD]: 

n^K^ l + 0.2271Afeff' ^ ' 

Because of the second effect, WMAP5 data favors A"cff = 3.046 over A'off = at the 
99.5% confidence level 11 . 

Neutrinos with mass ^ 1 eV become non-relativistic after the epoch of 
recombination probed by the CMB, so that allowing massive neutrinos alters matter- 
radiation equality for fixed ilmh^- Their radiation-like behavior at early times changes 
the expansion rate, shifting the peak positions, but this is somewhat degenerate with 
other cosmological parameters. Therefore, WMAP5 alone constrains "^rrii, < 1.3 
eV at the 95% confidence in a fiat ACDM universe, and this constraint relaxes to 
^ruy < 1.5 eV for a fiat wCDM universe, assuming Nefj = 3.046 [11]. Figure 17 of 
PI shows that there remain degeneracies between ^ and both Hq and erg ; in this 
paper we show that currently available low-redshift measurements can break these 
degeneracies, yielding improved and robust constraints on ^m^. 

After the neutrinos become non-relativistic, their free-streaming damps power 
on small scales and therefore modifies the matter power spectrum in the low redshift 
universe. The latest large-scale structure constraints from the Sloan Digital Sky Survey 
Luminous Red Galaxy sample in combination with WMAP5 yield '^m„ < 0.62 eV in 
a fiat ACDM model with Acff — 3.046 [12 . In addition, because the measurement of 
the matter power spectrum provides a separate constraint on the transfer function, this 
data combination also yields a constraint on A"cff = l.Sj^} 7 ^ ACDM cosmology 
with no massive neutrinos. We will show that the combination of low redshift data 
used in this paper provides tighter neutrino constraints in both of these cases. 
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Relaxing the assumptions about the cosmic expansion history, e.g., requiring 
w = —1 vs. allowing w as a free parameter [21 113) or allowing coupling in the 
dark sector can relax constraints on neutrino masses [5]. As an example we study 
constraints on neutrino masses in the model of |14| . allowing fife, w, and coupling 
strength ^ to vary simultaneously. 

In Section [2] we discuss the low-redshift observational constraints we will make 
use of to improve neutrino constraints, paying careful attention to how the constraints 
generalize to the models we consider. In Sections [3] and |4] we present new constraints 
on neutrino masses and the effective number of relativistic species at matter-radiation 
equality. We also make a detailed comparison with results from different cosmological 
probes available in the literature. In Section [5] we summarize our findings and argue 
that they are robust to a wide variety of alterations of the cosmological model. 

2. DATA AND METHODS 

In subsections 12.11 and 12.21 we discuss the low redshift constraints we make use of in 
our analysis, paying particular attention to their generalization to models not initially 
considered in their analyses. In subsection 12.31 we review our standard Bayesian 
methodology for computing parameter constraints and upper limits, and introduce 
the profile likelihood as a means to check the dependence on these constraints on our 
priors. 

2.1. Ho 

Ref. [TS] presents a redetermination of the Hubble constant that makes use of Cepheid 
variables in SN la host galaxies and in the "maser galaxy" NGC 4258, along with 
an expanded SN la sample at z < 0.1. We reinterpret the reported constraint, 
Hq — 74.2±3.6 km s^^ Mpc~^, as a constraint on the inverse luminosity distance at the 
effective redshift z = 0.04 (Riess, private communication). This correctly accounts for 
the slight dependence on w shown in Figure 14 of [15, , and generalizes the constraint 
to an arbitrary cosmology. In practice, when combining with other data sets, the 
cosmology dependence of this constraint is quite small. 

2.2. maxBCG 

Ref. [H] derives cosmological constraints from the maxBCG cluster catalog [T7], which 
identifies spatial overdensities of bright red galaxies in SDSS DR4-1- imaging data and 
results in a nearly volume-limited catalog with photometric redshifts in the range 
z = 0.1 to z = 0.3. MaxBCG clusters are classified by A^2007 the number of red- 
sequence galaxies within a scaled radius such that the average galaxy overdensity 
interior to that radius is 200 times the mean galaxy density. The group multiplicity as 
a function of A*2oo is one of two primary observable inputs to the maxBCG posterior. 
The second is the mean mass measured via gravitational lensing for 5 bins in iV200 
from [IH]. 

The analysis of Ref. [16] assumes a flat ACDM cosmology with N^ff — 3.046 
massless neutrinos, where ag and Q,m are allowed to vary, and h — 0.7, Ug — 0.96, 
and Q.bh^ = 0.02273 are held fixed. There are four nuisance parameters: two to define 
the mean power law relationship between N2Q0 and halo mass M, one to quantify the 
size of the log-normal scatter about this relation, and one for the lensing mass bias 
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parameter. The principal cosmological constraint derived from the maxBCG catalog 
is 

erg (an/0. 25)° ''i ^ 0.832 ± 0.033. (2) 

In Section 5.2.1, Ref. [16] addresses the limitations of their assumptions on 
cosmological parameters. They find that using a Gaussian prior on h — 0.7 ± 0.1 
and Us = 0.96 ± 0.05 does not alter their (Tg^^m'*^ constraint. In addition, for a prior 
with < 1 the agV,^^ constraint is not altered. Massive neutrinos suppress 

the linear power spectrum, but because as corresponds to the mass scale to which 
maxBCG is most sensitive, the impact of the alteration in power spectrum shape is 
minimal on the final constraints. In principle, massive neutrinos also alter growth 
in the nonlinear regime. Ref. [19] shows that power suppression is enhanced in the 
nonlinear regime for massive neutrinos, implying that the upper bound reported in the 
present paper is robust. Therefore, since the range of total neutrino mass we explore 
in this paper is well below 1 eV, we can safely apply the maxBCG constraint. In the 
case of massive neutrinos, we also study the maxBCG constraint in combination with 
WMAP5 for a ACDM model including spectral index running. When we evaluate 
the spectral index ns{k) for all k < 0.4 h/Mpc from the WMAP5 constraints in this 
model, the value falls between 0.91 and 1 for 68% of the time. This range satisfies the 
stated maxBCG tolerance to variation in the spectral index, and therefore we may 
apply the maxBCG constraint unaltered in the models we study including massive 
neutrinos. 

The maxBCG analysis did not directly address models in which iVoff is varied. 
However, we have verified that within the space of parameters allowed by WMAP5 
and the Hq constraint of 15 in a ACDM+A'cff model, the impact of A'off on the 
shape of the linear power spectrum can be compensated by a change in rig that is 
within the maxBCG stated tolerances. Therefore we may apply the constraint on the 
power spectrum amplitude through erg in Equation [2] without modification when A'cfi 
is varied. 

2.3. Markov Chain Monte Carlo and Profile Likelihoods 

We use the standard Markov Chain Monte Carlo (MCMC) scheme to probe the 
posterior distribution for several different cosmological models [11] [20l [21]. The 
publicly available WMAP5 chains which we resample have uniform priors on the 
parameters {i}bh^,^ch^,flA,T,ns,A'^,Asz}, the baryon density, cold dark matter 
density, dark energy density, reionization optical depth, scalar spectral index at 
ko = 0.002 Mpc~^, amplitude of curvature perturbations at fco, and SZ marginalization 
factor. The MCMC chain we produce for the the dark coupling model of [Tl] uses the 
same parameters, but replaces Vl^ with the parameter 0, the ratio of the approximate 
sound horizon to the angular diameter distance, which is used by CosmoMC [H] . 

For a more "frequentist" statistic, the profile likelihood can be computed from 
the MCMC chain [23]. Since we have 7 or more uninteresting parameters and only 
one parameter for which we want to report constraints (say, /?), for each value of 
/3 we find the maximum likelihood value Lp^ regardless of the values assumed by 
other parameters. We then consider hiLp / Lmax as a function of the model parameter 
^, or function thereof; here Lmax is the maximum likelihood over the entire chain. 
We examine the profile likelihood for both ^ (Section 13.31) and N^s (Section 
14. 2p . As long as the posterior distribution is singly peaked, this profile should be 
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Figure 1. ACDM WMAP5-only constraints with "^rrii, = are shown in blue. 
The black transparent contours show how the constraints on fTg(f2m/0.25)'''*^ 
and Ho degrade when ^m,j is left as a free parameter. The dotted lines show 
contours on which ( m^) = 0.2, 0.4, 0.6, ...,1.4 eV for the WMAP5-only posterior 
distribution. The green shows the constraints when maxBCG and Hg constraints 
are also included. The 95% confidence upper limit is reduced from 1.3 eV to 0.4 
eV. 



independent of the Bayesian priors; but note that the profile hkelihood will become 
noisier as InL/L^ax decreases, since the MCMC method means the chains do not 
densely sample low \a.L/ L^ax regions. 

3. CONSTRAINTS ON THE SUM OF NEUTRINO MASSES 

We begin in Section 13.11 by studying the ACDM model and a few one-parameter 
extensions explored in the WMAP5 analysis. In section 13.21 we study the dark 
coupling model of 14 as an example of a model which significantly relaxes the model 
assumptions about the expansion history by varying three additional parameters. We 
consider the profile likelihood to study the dependence of our constraints on our 
cosmological priors and check the consistency of the data sets we examine. Finally we 
compare with previous constraints on ^ to^ in the literature using a variety of data 
sets. 

3.1. ACDM and one-parameter extensions 

Figure [T] shows the WMAP5-only constraints in a ACDM model when ^ = 
(blue contours). The solid black contours show that when ^ rui, is allowed to vary, the 
constraints on Ho and cr8(i7,„/0.25)°'^^ weaken substantially. While both the maxBCG 
and Hq constraints are in good agreement with the ACDM values, they effectively 
exclude high neutrino masses (green contours). Adding both constraints reduces the 
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95% confidence limit on '}2,m,^ from 1.34 to 0.4 in a flat ACDM cosmology. Figure [5] 
shows the same scenario, but when BAO and SN data have also been combined with 
WMAP5. Since those probes break the degeneracy between Hq and in ACDM 
cosmology, Hq is well constrained, almost independently of '^m^. In this case the 
maxBCG result drives the reduction of the 95% confidence bound from 0.67 eV to 
0.31 eV. 

In Table [T] we present these and several other one dimensional constraints on 
^ my for the ACDM model, as well as several other data and model variants. In 
the upper half of the table the limits are Bayesian, i.e. they represent the ^ nii, 
value below which 95% of the weighted points in the MCMC chain lay. The lower 
half of the table approximates a frequentist interpretation (a more detailed discussion 
of this is in sec. 13.31) : we report the largest value of '^m^, in the chain for which 
\nLj2m,jLmax > "2. While this limit can be interpreted in terms of numbers of 
standard deviations only for Gaussian distributions, it is prior-independent. We have 
checked that the ^ rrii, distribution is well approximated by a Gaussian. 

The constraint on ^ from WMAP5-|-maxBCG-|-7fo does not relax when 
running or tensors are allowed. The constraint relaxes slightly to '^rui, < 0.47 eV 
at the 95% confidence level when w is allowed to vary, but BAO and SN constraints 
are also included. The constraint relaxes because the maxBCG constraint - J2 "^i' 
degeneracy is less steep for this model, since changing w also impacts the growth 
of structure. These examples indicate that when the expansion history is fixed, 
^ < 0.4 eV (95% confidence) is a robust upper limit; changes in the primordial 
fluctuations that are consistent with WMAP5 do not impact this bound. However, 
when the expansion history is altered at late times (w ^ —1 but constant), this bound 
relaxes slightly. As an extreme example of this, we study the dark coupling model of 
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Bayesian and Frequentist upper 95% C.L. bound on 



model 


base dataset 


- 


+maxBCG 


+Ho 


+niaxBCG+ffo 


ACDM 


WMAP5 


1.3 


1.1 


0.59 


0.40 


ACDM 


WMAP5+BA0+SN 


0.67 


0.35 


0.59 


0.31 


ACDM +a 


WMAP5 


1.34 


1.25 


0.54 


0.39 


ACDM +r 


WMAP5 


1.36 


1.18 


0.83 


0.40 


wCDM 


WMAP5+BA0+SN 


0.80 


0.52 


0.72 


0.47 


dark coupling 


WMAP5+A,,o(fc)+SN 






0.51 




ACDM 


WMAP5 


1.3 


0.82 


0.50 


0.33 


ACDM 


WMAP5+BA0+SN 


0.66 


0.32 


0.56 


0.30 


ACDM +a 


WMAP5 


1.28 


1.15 


0.63 


0.43 


ACDM +r 


WMAP5 


1.23 


0.86 


0.72 


0.30 


wCBM 


WMAP5+BA0+SN 


0.80 


0.44 


0.74 


0.44 


dark coupling 


WMAP5+Aa;o(fc)+SN 






0.55 





Table 1. 95% upper limit on the sum of the neutrino masses "^mi, in eV. In 
the upper half of the table the limits are Bayesian, i.e. they represent the nii, 
value below which 95% of the weighted points in the MCMC chain lay. The lower 
half of the table approximates a frequentist interpretation: we report the largest 
value of ™ the chain for which In L-^^^^ / Lmax > —2. For the first five 

models we use the WMAP5 public MCMC chains to compute the upper bounds; 
these chains use WMAP5 year data, and where noted, BAO constraints from I24 | 
and the Union supernova sample 25 . Column 1 specifies the cosmological model 
assumed. All five of these models assume flatness, and w is assumed constant if 
varied, a = dris/dlnk indicates running in the primordial power spectrum and r 
indicates allowing for primordial tensor modes. Column 2 specifies the original 
datasets used to determine the constraints shown in Column 3, and Columns 4-6 
give the resulting upper bounds after applying additional observational constraints 
from maxBCG (Column 4), Hq, (Column 5), and both of these (Column 6). In the 
last column we include the constraint from Section 13.21 on the model of Ref . I14| . 
In that model, curvature, w, and dark coupling ^ are simultaneously varied. 



[14] next. 

3.2. Dark coupling and massive neutrinos 

In this subsection we consider the particular model of coupling in the dark sector 
proposed by Ref. [13] . In their model the densities of dark matter and of dark energy do 
not evolve independently as the universe expands but are coupled, with an interaction 
strength parameterized by ^ as follows: 

V.r('^,„).=^i?P.e«r (3) 

V^r^'^rfe). = - ^Hp,,ut"\ (4) 
where T/^ , and T/^, , are the energy momentum tensors for the dark matter and 

dark energy components, respectively, u"^"^ denotes the dark matter four velocity and 
H the Hubble parameter. The coupling strength ^ alters the expansion history through 
changes in the matter and dark energy density evolution [14] : 



Pd,„(a) - pl^a + p^e^(l - a""" )a-' (5) 

p,M) =pL«-''^+'""^ (6) 



Robust Neutrino Constraints from Low Redshift Observations 



8 



where w'^^ = w + ^/3, the superscript denotes the present-day quantity, and a 
denotes the scale factor. We adopt this model as a suitable example for which to study 
whether the constraint on ^ relaxes when several free parameters are allowed to 
simultaneously vary. It was shown by Ref. [IHl 123 dH] that allowing for a coupling 
in the dark sector dramatically relaxes cosmological constraints on neutrino mass. 
Ref. |14j presents a slightly different model for dark coupling, which improves on 
existing models in the literature as it does not suffer from instabilities; they present 
cosmological constraints on this dark coupling parameterization in a very general 
model, allowing the curvature, equation of state of dark energy w > —1, total neutrino 
mass, and the amplitude of dark coupling ^ < between dark matter and dark energy 
to vary along with the standard ACDM cosmological parameters. These additional 
parameters allow the expansion history and growth of structure to deviate appreciably 
from the ACDM assumption. We will explore the cosmological constraints on the 
coupling strength ^ in a future work alongside this model's non-standard growth of 
structure; because of the large deviations from the ACDM growth of structure, the 
maxBCG cluster constraint, which is based on the standard spherical collapse-based 
correspondence between the linear matter power spectrum and the resulting halo mass 
function in ACDM cosmologies, cannot be applied to the [Ml model. For this model, 
we therefore substitute the maxBCG constraints with the SDSS DR7 Luminous Red 
Galaxies Phaio{k) constraint of [H] as a robust low redshift probe of X)"^!" 

We study the dark coupling model using the MCMC method as implemented 
in the COSMOMC package [22]. Several datasets must be included to avoid long 
degeneracies in the parameter space. We include WMAP5 data along with several 
low redshift datasets: the new Hg constraint [15], SDSS DR7 Luminous Red Galaxies 
halo power spectrum Phaioik) |12| . and Union SN sample [25j . The only important 
prior we enforce is ftdmh^ > 0.001; surprisingly, even models consistent with these four 
cosmological probes can have negligible fldmh^ today since the coupling accelerates 
the transition to dark energy domination. Note that massive neutrinos suppress 
Phaioik) in a scale-dependent fashion, thus, by including the Phaioik) measurements, 
one expects to be able to disentangle the effects of neutrino mass from the effects of 
alterations of the late-time expansion history. Figure [3| shows constraints from these 
four observational probes (red contours): is not degenerate with dark coupling 

^, and the resulting one dimensional constraint on ^ m^, is 0.51 eV. For completeness, 
we also show in Figure [3]how the constraints would tighten if the maxBCG constraint 
were applicable (blue contours). The constraint on the neutrino mass is unaffected. 

3. 3. Comparison with constraints from the profile likelihood 

We first verify that all three primary datasets are consistent with each other in 
the ACDM cosmology by comparing the maximum likelihood value of the WMAP5 
data alone with the maximum likelihood after adding the additional observational 
constraints. We find that \nLmax decreases by 0.5 with the addition of the two 
independent Gaussian priors from the Hq and maxBCG constraints; the decrease in 
In Lmax is the same when we add massive neutrinos to the model. Therefore, our tight 
constraints are not arising from tension between the datasets we are using. To compute 
the profile likelihood as described in Section[231 we find the maximum likelihood value 
over the MCMC chain in bins of X] '^ii/. Figure |4| shows that this function is well- fit by 
a Gaussian in each case. We can therefore tentatively interpret the largest ^ m^, value 
for which In L/Lmax > ^2 as a prior- independent ~ 2a upper bound on '^m^. We 
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Figure 3. Two-dimensional constraints in ^ — '^rriy in the dark coupling 
model of ref. I14| . The red contours are the results from the CMB-I- Hq 1151 + 
PhaloiJ^) 1121 + Union SN |25| . With this data-set combination there is no longer a 
degeneracy between coupling and neutrino mass. As discussed in Section 13.21 the 
strong changes in the growth of structure caused by the dark coupling prevent us 
from applying the maxBCG constraint to this model. However, the blue contours, 
which also include the maxBCG constraint, demonstrate that the constraint on 
is not tightened even if this additional constraint could be included. 

report this value in the bottom half of Table [TJ in general we find excellent agreement 
with the Bayesian results. 

3-4-. Comparison with other constraints 

In Table [2] we compile other ^ mi, constraints available in the literature and compare 
with the results we have found. WMAP5 provides a robust upper limit in a ACDM 
model, < 1-3 eV; allowing w ^ — 1 relaxes the constraint only slightly to 

Yyi^ty < 1.5 eV. In a study of WMAP3 data, Ref. [29] find that correcting the 
likelihood approximation on large angular scales as well as accounting for an over- 
subtraction of point sources affecting small angular scales, the WMAP3-only 95% 
upper bound improved from 1.9 to 1.57 eV. '^m^, remains degenerate with both 
and (78 within the WMAP5 constraint. Using BAO and SN reduces the degeneracy 
with Hq, reducing the limit to 0.67 eV (w — —1) or 0.8 eV when w is also allowed to 
vary. Using both the /c-dependent suppression of power and BAO feature in the matter 
power spectrum derived from the SDSS Luminous Red Galaxy sample, Ref. ^12j finds 
^TTii, < 0.62 eV when combined with WMAP5 for w = — 1. That analysis improves 
the nonlinear modeling by removing the additional shot noise contribution from 
galaxies occupying the same dark matter halo and marginalizing over the remaining 
modeling uncertainties. Ref. [30] compares the impact of two nonlinear modeling 
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0.0 0.5 1.0 1.5 2.0 

2 m„ (eV) 

Figure 4. The profile likelihood defined in Section l2.3l in bins of A( m,,) = 0.2 
eV for the ACDM model with WMAP5 only (crosses), WMAP5+maxBCG 
(stars), WMAP5+Ho (triangles), and WMAP5+maxBCG+Ho (diamonds). The 
black curves overlay a quadratic fit to these points, illustrating that a Gaussian 
curve provides a good fit to this one-dimensional distribution. 

approaches on the cosmological inferences for several extensions of standard ACDM 
cosmology, and advocates marginalizing over an additional shot noise term. 

Additional progress requires probing the amplitude of mass fluctuations. This 
has been done with several samples probing the cluster mass function. Ref. |31] 
derives a low and high redshift cluster mass function from Chandra observations of 
clusters in two ROSAT surveys. Their measurement provides a constraint from the 
normalization of the local mass function o-8(f^m/0.25)° '^'^ = 0.813 ± 0.013 ± 0.024 
(statistical and systematic errors, respectively); this is similar to and in agreement 
with the maxBCG constraint we have used in our analysis. Additional cosmological 
information is available by comparing the high and low redshift mass functions. Their 
tighter statistical constraint on the local mass function normalization translates to 
tighter constraints on ^ to^ < 0.33 eV, but with the caveat that ^ 0.1 eV of systematic 
error has not been included. Moreover, Ref. [32] points out that the analysis of [3T] 
does not fit for both the cluster scaling relations and cosmological parameters self- 
consistently. 

Ref. [33] uses the cluster mass function from weak lensing to derive a similar 
cluster normalization constraint: cr8(ilm/0.25)" '^^ — 0.72j^Q Qg. Because of both the 
larger errors, lower central value, and combination with WMAP3 rather than WMAP5, 
they find a relatively weak constraint, "^m^ < 1.43 eV. 

Because so far only upper limits can be obtained from cosmology on the neutrino 
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masses, the tightness of the ^ constraints depends drasticahy on where the lower- 
redshift trg determination intersects the CMB ^ nii, — erg degeneracy. An example of 
this sensitivity is the constraint obtained including the Ly-a forest power spectrum 
(LYA) measurements |34j . In that work, ~ 2a tension between the power spectrum 
amplitudes preferred by WMAP3 and LYA when ^ nii, = was reported, with LYA 
preferring a larger amplitude. This translates into an extremely tight constraint 
on ^ mi, < 0.17 eV when LYA is combined with many other observational probes 
(WMAP3, SN, SDSS BAO, SDSS and 2dF P{k), and small scale CMB experiments). 
With updated datasets for WMAP and SN but including fewer datasets altogether, 
and in addition allowing for running of the spectral index, we find that LYA data still 
reduces the constraint to X]"^^ < 0-28 eV. Moreover, Ref. [35] finds possible evidence 
for an inverted temperature-density relation, which would change the interpretation of 
the LYA power spectrum measurement and significantly weaken the Ref. [34j bound. 

Table [1] demonstrated the improvement on ^ m,^ when tight Hubble constant 
measurements are included [TS]. Refs. use the stars in passively-evolving 

red galaxies as cosmic chronometers to obtain a cosmological-model independent 
measurement of H{z) from z = to z 1; adding these additional constraints on 
H{z) also reduces the constraint on ^ rrii, to 0.5 eV. Weak lensing provides additional 
constraints on the amplitude of fiuctuations on small scales which can be used to 
constrain "^my; both Ref. [55] and Ref. [SS] find a limit of ~ 0.54 eV. Finally, there 
have also been constraints derived from luminosity-dependent biasing [JO] [41] . The 
most recent work by [41] combines WMAP5 with the [42] LRG P{k) and luminosity- 
dependent bias measurements to find '^rriij < 0.28 eV in ACDM; this constraint 
relaxes to 0.59 eV in a ?z;CDM cosmology. While the modeling procedure can be quite 
complex, this approach is certainly competitive with other probes. 

In conclusion, WMAP5 provides a robust 95% confidence upper bound of ~ 1.5 
eV. Probes of the redshift-distance or H{z) relation break key degeneracies and reduce 
the constraints to ^ 0.5 — 0.8 eV. Similarly, adding the measurement of the galaxy 
power spectrum shape for SDSS DR7 [12] reduces the bound to 0.62 eV. Beyond this, 
probes of the amplitude of clustering as through measurements of the cluster mass 
function, weak lensing, LYA, or luminosity-dependent clustering are necessary. The 
most optimistic interpretation of the currently available data reduces the bound to 
~ 0.3 eV. All of these methods require more complicated modeling and need further 
quantification of the remaining systematics. 

From current solar and atmospheric constraints on neutrino mass differences 
Ato^, three types of neutrino mass spectra are allowed: the normal hierarchy 
(mi ^ m2 <C ma), the inverted hierarchy (m^ <C mi ~ m2) and the quasi-degenerate 
hierarchy (mi m2 ^ ms). Constraints on the absolute neutrino mass scale would 
greatly help in distinguishing there three cases as indicated in Figure|5j Direct neutrino 
mass searches based on tritium beta decay so far have imposed an upper limit of 2.3 
eV at the 95% CL but to distinguish e.g. the quasi-degenerate hierarchy from the 
other cases a mass sensitivity of ~ 0.2 eV is needed. The neutrino mass hierarchy is 
one of the unknowns within the neutrino sector and its extraction would be crucial for 
determining the neutrino character (Dirac or Majorana) and measuring leptonic CP 
violation. A determination of the neutrino mass hierarchy would also rule out some 
neutrino mass models based on lepton flavor symmetries and/or grand unification 
schemes [IJ . If future cosmological data finds no evidence for a neutrino mass down 
to 0.1 eV, that would imply that the neutrino mass hierarchy is normal, allowing, 
in principle, a clean measurement of CP violation in the leptonic sector (provided 



Robust Neutrino Constraints from Low Redshift Observations 



12 




Figure 5. vs. the mass of the lightest neutrino in the normal and inverted 

hierarchies. 



the mixing angle ^13 is not very small), see Ref. |45]. If future cosmological data 
finds a positive signal above 0.1 eV it would be impossible to measure unambiguously 
the hierarchy with cosmological information only. However, if such a signal is found, 
direct neutrino mass searches, as tritium beta decay experiments, could infer the 
mass of the lightest neutrino and provide the key to determine the hierarchy and 
leptonic CP violation. Moreover, if the neutrino mass hierarchy turns out to be 
inverted, next generation neutrino-less double beta decay experiments could extract 
the neutrino character (for a review of the expected sensitivities, see |46]). A Majorana 
neutrino character, together with a non-vanishing leptonic CP violation would point 
to leptogenesis as the responsible mechanism for the matter-antimatter asymmetry of 
the universe. 

4. Constraints on the number of relativistic species A'rci 

We now consider constraints on the effective number of neutrino species: in this case 
we assume that neutrino masses are negligible and in the analysis we keep ^ m^, fixed 
to 0. While in principle both N^-d and ^ m,^ should be varied simultaneously, in 
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95% upper limits on ^ m,y from the literature 



model 


base dataset 


U.L. (eV) 


sys. errors 


Ref. 


ACDM 


WMAP5 


1.3 


m 


m 


wCBM 


WMAP5 


1.5 


m 
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ACDM 


WMAP5+BA0+SN 


0.67 


- 




wCDM 


WMAP5+BA0+SN 


0.80 


- 




ACDM 


WMAP5+F,,a/o(fc) 


0.62 


marginalized 


112 


wCDM 


WMAP5+BA0+SN+ 










Chandra Clusters 


0.33 


±0.1 


m 


ACDM 


WMAP3+WL mass function 


1.43 




m 


ACDM 


WMAP3+Ly-a+CMB small 










+BAO+SN+SDSS & 2dF P{k) 


0.17 


IGM [35] 


m 


ACDM 


WMAP5+Ly-a+CMB small 








+a 


+Ho 2001 47 


0.28 


IGM [35] 




ACDM 


WMAP3+Ho+Hiz) 


0.50 






ACDM 


WMAP5+SN+BA0+WL 


0.54 


+0.04 


[38 , [39 


ACDM 


WMAP5+b{L)+PLRG{k) 2006 


0.28 


bias modeling 


m 


wCBM 


WMAP5+6(L)+PL7?G(fc) 2006 


0.59 


bias modeling 


m 



Table 2. 95% upper limit on the sum of the neutrino masses rrii, in eV 
reported in other studies. All constraints include WMAP, and have been sorted 
by horizontal lines by the type of secondary probe: geometrical probe(s) BAO 
and/or SN, galaxy clustering, cluster mass function, Lyman-« forest, cosmic 
chronometers, weak lensing, and luminosity-dependent galaxy biasing. In the 
second Lyman-a analysis listed, a = dUg/dlnk is also a free parameter. 



practice the two are no longer degenerate, and so the assumption ^ = docs not 
change the constraints much [48 . 

Equation [1] defines the major degeneracy between A^id and ^Imh^ in the WMAP5 
data. With the strict assumptions of a ACDM cosmology, the measured angular 
diameter distance can only be recovered if Hq varies with N^ci, while the rim and 
fifcft,^ constraints do not relax when TVrei is allowed to vary. This explains the strong 
degeneracy between iVrei and Hq in the WMAP5 posterior distribution. There is also 
a weaker degeneracy between iVrei and rig. While WMAP5 tightly constrains the 
amplitude of fluctuations on very large scales, fco ^ 0.002 Mpc~^, within the WMAP 5 
degeneracy and assuming a power law primordial spectrum, the allowed linear P{k) 
pivots about fco as Nj^^i varies, so that the amplitude of power near the peak at fc ~ 0.03 
Mpc~^ (and beyond) increases with iVj-ei. This tilt causes the degeneracy between A^'rei 
and CTg , and can be constrained both by measuring the galaxy power spectrum shape 
and constraints on the cluster mass function. Finally, since flm and fli,h^ are both 
tightly constrained, the variation of Hq with iVici induces a degeneracy between TViei 
and rih/flm- This could be probed through gas fractions in clusters or eventually 
through the amplitude of the BAO signal in large scale structure data. These well- 
defined one-dimensional degeneracies would degrade in a less-restricted model of the 
expansion history, so that further geometrical constraints would be needed to constrain 
Nj-ei and other parameters like flk and w simultaneously. 
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68% and 95% confidence intervals for N^d 
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•5- ' ' -0.67V-1.24J 


WMAP5 








4.o(i}j) 


3.76(11;}?) 


WMAP5+BA0 
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4.13(11:^7^) 


3.73(t}j^) 


WMAP5 

+ Phalo{.k) 








4.i6e}i?) 
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Table 3. The upper half of the table shows the Bayesian mean, 68% and 95% 
confidence intervals on N^-^i for a ACDM cosmology. Here HST refers to the 
I47| constraint on Hq : 72 ± 8 km s~^ Mpc~^. When the |15l Ho constraint is 
included in columns four and five, the 1471 HST prior is undone. As in Table [T] 
in the lower half we compute the maximum JV^ei for which InL/Lmax > —2; for 
ease of comparison we reproduce the central values from the upper portion of the 
table. Since the distributions are significantly non-Gaussian unless WMAP, Hq, 
and maxBCG constraints are included, these "frequentist" constraints cannot be 
interpreted as confidence intervals. 



4-1- Nrci in ACDM using several datasets 

Figure [6] shows the degeneracies in the WMAP5-|-P;ia/o(fc) posterior distribution 
between A'roi and both Hq and the maxBCG-constrained parameter cr8(51m/0.25)° '*^; 
the results are qualitatively the same when the Phaioik) constraints are excluded. 
Table [3] illustrates constraints on TV^ei in ACDM for three different base datasets: 
WMAP5 alone, WMAPS-t-BAO-hSN-hg^ HST prior {Hq = 72 ± 8 km s^^ Mpc"!), 
and WMAP5-|-P;iaio(fc) from [T^]. Both the maxBCG and updated Hq measurement 
reduce the error on iVroi significantly in all three cases. When both of these constraints 
are taken together, additional data sets other than WMAP5 do not shift the mean 
value of iVrci or reduce its error: N-^ci = 3.76lg'gg. Therefore the maxBCG and 
Hq constraints are sufficient to break the key degeneracies in the WMAP5 data in 
this model. Since the geometrical constraints from BAO and SN do not contribute 
significantly to the N^c\ constraint, they would probably be useful in constraining more 
general models (i.e., models which vary curvature or the behavior of dark energy) in 
the context of freely varying A^roi- 

4-2. Profile likelihood for A^roi *^ ACDM model 

Figure m and Table [3] indicate that the likelihood surface is skewed - contours are more 
widely separated at large N^ci- This finding is in line with the Ref. [23] result that 
the posterior of iVrci only approaches a Gaussian distribution when several datasets 
are combined. In Figure [7] this skewness is evident in the profile likelihood for Nrei- 
A Gaussian provides a reasonable fit to the profile likelihood only when both Hq and 
maxBCG constraints are included. Therefore, the "frequentist" constraints reported 
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Figure 6. ACDM WMAPb+Phaioik) constraints on A^rei vs. Hq and tiie 
maxBCG constraint, (Tg(r2m/0.25)'^'*"'^ (black), and after importance sampling 
the WMAP5+Pjia;o(A:) constraints using the |15| constraint on Ho (blue), the 1161 
maxBCG cluster constraint on (Tg(f!m/0.25)'' '*^ (red), and both (green). These 
contours arc similar when the Phalo{k) constraints are excluded. 
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Figure 7. The profile likelihood defined in Section 12.31 in bins of 
AiVrci = 0.4 for the ACDM model with WMAP5+Ph„io(fc) only (crosses), 
WMAP5+Phaio(fc)+maxBCG (stars), WMAP5+Phaio{k)+Ho (triangles), and 
WMAP5+Phaio(''^)+™a'XBCG+Ho (diamonds). The curves are similar when 
Phalo{k) is excluded but Hq and/or maxBCG constraints are included. The 
distributions are significantly non-Gaussian, though when both the maxBCG 
and Hq constraints are included, a Gaussian distribution (solid curve) provides a 
reasonable fit to the profile likelihood. 

in the lower half of Table [3] are only meaningful for the last column, which includes Hq 
and maxBCG constraints. In that case, there is good agreement with the Bayesian 
95% confidence interval. 

4-3. Comparison with other N^-d constraints 

Table H] shows constraints on N^^ from other data in the literature. While WMAP5 
provides robust evidence for a lower bound on iVioi, by itself it cannot provide an 
upper limit. The strong correlation between iVj-ei and Hq in the ACDM model 
makes direct probes of Hq and H{z) the most effective way to constrain iVroi with 
low redshift observations. The new HST constraint clearly provides a significant 
improvement over WMAP5 in combination with previously available geometrical 
constraints (WMAP5+BA0+SN+HST). The additional constraints on the expansion 
history from [36 improve even further, with the WMAP5+i?o error on A'roi of 1.0 
being reduced to 0.5 with the combination WMAP5+Ho+H{z). The Ly-a constraints 
are also quite promising as a probe of the power spectrum amplitude at large fc, but 
should be reevaluated with the new data and in the context of the possible systematics. 
Comparison of measured isotope abundances with BBN calculations provide still the 
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base dataset 




systematic errors 
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Table 4. Constraints on A^j-el reported in other studies. Here HST refers to the 
[471 constraint on iJo : 72 ± 8 km s"! Mpc-^. All constraints include WMAP 
except the Big Bang Nucleosynthesis (BBN) constraints, and have been sorted 
by horizontal lines by the type of secondary probe; geometrical probe(s) BAO 
and/or SN, galaxy clustering, primordial abundances from BBN, Lyman-a forest, 
and cosmic chronometers. 



tightest and completely independent constraints on A^ici. 
5. Conclusions 

In this paper we have demonstrated that probes of the cluster mass function and 
increased precision on the Hubble constant can break key degeneracies with CMB 
observations and yield excellent constraints on both the number of species and sum of 
the masses of cosmological neutrinos. When the expansion history is fixed to ACDM, 
current constraints on Hq ^15] and the cluster mass function [16] constrain ^ m„ < 0.4 
eV at 95% confidence. This bound relaxes to 0.5 eV in the two extended models we 
considered: wCDM, and the dark coupling model of [2], which allowed curvature, 
w, and coupling strength ^ to vary. Probing the mass function using X-ray clusters, 
Ref. [31] combines X-ray cluster data with WMAP5-hBA0-hSN and finds ^2^,. < 0.33 
eV in a wCDM cosmology, though the systematic errors still must be clearly quantified. 
The optimistic upper bound available from current data, '^m^, ~ 0.3 eV, almost 
excludes the scenario in which the neutrino masses are quasi-degenerate. 

The constraint on the number of relativistic species iVioi from probes of the cluster 
mass function has not been widely explored so far. We find that the combination of 
WMAP5, [H] Hq, and the maxBCG cluster mass function constraint provides an 
excellent constraint: iVi-oi = 3.76lQ'gg. However, we point out that this constraint 
does not improve when BAO and SN data are also included; those data sets have 
been shown to have excellent constraining power on both and w [H]. Therefore, 
we may expect that the constraint on N,-ci will not relax dramatically if we allow for 
these additional parameters but include all of these datasets, particularly because of 
the complementarity of the growth of structure and expansion history constraints we 
have employed. These constraints are not competitive with the most optimistic errors 
from BBN [7], but provide an important consistency check using probes significantly 
after the BBN epoch. 
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